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ABSTRACT 

We describe the motivation, design and implementation of the CORNISH survey, an arcsecond resolution 
radio continuum survey of the inner Galactic plane at 5 GHz using the Karl G. Jansky Very Large Array 
(VLA). It is a blind survey co-ordinated with the northern Spitzer GLIMPSE I region covering 10°< I <65° 
and |6| <1° at similar resolution. We discuss in detail the strategy that we employed to control the shape of 
the synthesised beam across this survey that covers a wide range of fairly low declinations. Two snapshots 
separated by 4 hours in hour angle kept the beam elongation to less that 1.5 over 75% of the survey area 
and less than 2 over 98 % of the survey. The prime scientific motivation is to provide an unbiased survey for 
ultra-compact H ii regions to study this key phase in massive star formation. A sensitivity around 2 mJy will 
allow the automatic distinction between radio loud and quiet mid-IR sources found in the Spitzer surveys. 
This survey has many legacy applications beyond star formation including evolved stars, active stars and 
binaries, and extragalactic sources. The CORNISH survey for compact ionized sources complements other 
Galactic plane surveys that target diffuse and non-thermal sources as well as atomic and molecular phases 
to build up a complete picture of the ISM in the Galaxy. 

Subject headings: Astrophysical Data, Stars, ISM 
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1. Introduction 

The study of star formation, stellar populations and 
Galactic structure is currently being transformed by a 
new generation of Galactic plane surveys. These surveys 
have much higher resolution and sensitivity than previ- 
ous surveys and are at infrared and longer wavelengths 
that can penetrate the high extinction of the Galactic 
mid-plane. They also cover sufficiently wide areas of 
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the plane to provide samples that are large enough to 
be representative and statistically robust for rare and 
short-lived phases of evolution. Blind surveys that cover 
contiguous, wide areas also allow the study of phenom- 
ena and trends across a large range of scales. A wide 
wavelength coverage is required to distinguish and char- 
acterise the different populations found at both the early 
and later phases of stellar evolution and so it is also im- 
portant to ensure that the same wide areas are covered by 
all the relevant wavelengths. Here we describe the design 
and implementation of a new radio continuum survey of 
the Galactic plane and its role in the context of the other 
multi- wavelength surveys. 

1.1. A New Generation of Galactic Plane Sur- 
veys 

Advances in various wide-field technologies have led 
to a new generation of Galactic plane surveys. Those 
with relevance here, and to star formation in partic- 
ular, are listed in Table [T] Leading the way is the 
GLIMPSE I (Galactic Legacy Infrared Mid-Plane Survey 
Extraordinaire) survey (Benjamin et al. 2003; Church- 
well et al. 2009). This Spitzer IRAC Legacy Pro- 
gramme covered the inner galaxy (10° < I < 65° and 
-65° < I < -10°, |6| < 1°) at 3.6 Aim, 4.5 /Ltm, 5.8 /xm 
and 8.0 /im. GLIMPSE is two orders of magnitude more 
sensitive (limit at L-band (3.6 fim.) « 14 — 15 magnitude 
although usually confusion limited) and has ten times 
higher spatial resolution (1.4 — 1.9") than any previ- 
ous mid-IR survey. It has catalogued over 49 million 
sources in a wavelength regime which preferentially se- 
lects sources with hot circumstellar dust emission such 
as young and evolved stars (e.g., Robitaille et al. 2008). 
More embedded and cooler sources are the subject of the 
MIPSGAL survey with Spitzer (Carey et al. 2009) and 
the PACS element of the Hi-GAL survey with Herschel 
(Molinari et al. 2010). 

An even deeper probe of the general stellar popula- 
tion is becoming available via the deep near-IR Galactic 
Plane Survey (GPS) that is part of the UK IR Deep 
Sky Surveys (UKIDSS) programme (Lucas et al. 2008). 
This survey is about 3 magnitudes more sensitive and 
two to three times better resolution than the 2MASS 
all sky near-IR survey (Skrutskie et al. 2006) and will 
detect upwards of a billion sources. The addition of 
near-IR data, where photospheric and scattered light 
contributions dominate, allows a much better separa- 
tion and characterisation of the general stellar popu- 
lation, than from mid-IR colours alone. Near-IR sur- 
veys are also being used to map the extinction due to 
molecular clouds (e.g., Rowles & Froebrich 2009). The 
optical Ha survey IPHAS (Drew et al. 2005) is pro- 
viding both stellar characterisation, 3D extinction maps 
and tracing nebular emission in less obscured regions of 
the northern plane with the southern VPHAS-f survey 
(see www.vphasplus.org) scheduled to do the same in the 
south. 

To place the stellar populations in a Galactic con- 
text, commensurate studies of the various components 



of the ISM are required. For the study of star for- 
mation in particular, the molecular component is pro- 
vided in part by the BU-FCRAO i^cQi-O Galac- 
tic Ring Survey (Jackson et al. 2006), which covers 
over half of the northern GLIMPSE area with exten- 
sions over much of the northern mid-plane (Mottram 
& Brunt 2010). The mostly optically thin "cO 1-0 
data traces the distribution and dynamics of the cold 
molecular gas. The cool dust in these molecular clouds 
can also be mapped via its sub-millimetre continuum 
emission. This is currently being achieved with the AT- 
LASGAL survey (Schuller et al. 2009) at 870 /^m and 
BOLOCAM survey at llOO^im (Aguirre et al. 2010), 
whilst the forthcoming JCMT Galactic Plane survey 
(JPS, www.jach.hawaii.edu/JCMT/surveys/) will map 
the continuum at 450 and 850 /zm with higher sensitiv- 
ity and resolution. Together with the SPIRE element 
of the Hi-GAL survey with the Herschel satellite (Moli- 
nari et al. 2010) this will provide temperature and dust 
emissivity information across the GLIMPSE region. 

The more widely distributed atomic hydrogen compo- 
nent is covered by the International Galactic Plane Sur- 
vey at about 1' resolution (www.ras.ucalgary.ca/IGPS/). 
This is made up of the VLA Galactic Plane Survey 
(VGPS, Stil et al. 2006), Southern Galactic Plane Sur- 
vey (SGPS; McClure-Grifiiths et al. 2005) and Cana- 
dian Galactic Plane Survey (Taylor et al. 2003). A 
sub-section of the plane is covered by the much deeper 
GALFA survey of H i at 4' resolution using the Arecibo 
telescope (Peek et al. 2011). The Hi component traces 
where most of the mass of interstellar gas is located and 
its velocity structure. 21cm absorption lines can also 
help to solve the near/far distance ambiguity in relation 
to the kinematic distances derived towards the embed- 
ded Hii regions (Sewilo et al. 2004; Fish et al. 2003) 
and dense, molecular clouds (Jackson et al. 2002; Gib- 
son et al. 2005; Busfield et al. 2006; Roman-Duval et al. 
2009). This combination of atomic and molecular kine- 
matic data can therefore provide the 3D Galactic setting 
of the neutral gas. 

The one major component missing from the multi- 
wavelength surveys of the Galaxy is that of the ionized 
gas. No existing radio continuum survey has sufficient 
resolution, depth or coverage to provide the data that 
would complete the picture of the Galaxy. Ionized gas 
arises in nebulae and stellar winds around hot stars and is 
crucial to understanding key phases of the early and late 
evolution of high and intermediate mass stars. Dense, 
photo-ionized regions around young and evolved stars are 
often optically thick with a turnover frequency of around 
a few GHz, where the radio continuum spectral index a 
{Si, oc ly") changes from the optically thick value of -1-2.0 
to the optically thin slope of —0.1. Ionized stellar winds 
have a positive spectral index a ^ -1-0.6, and so it is 
also much more efficient to carry out systematic searches 
for these sources at high frequencies. Furthermore, the 
resolution should be comparable to the arcsecond reso- 
lution of the IR surveys to allow IR counterparts to be 
uniquely identified. It was to fill this need for a high fre- 
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quency, high resolution radio continuum survey covering 
the area of the Spitzer survey that the concept of the 
Co-Ordinated Radio 'N' Infrared Survey for High-mass 
star formation (CORNISlfl) survey was conceived. 

1.2. Previous Radio Surveys 

Many single-dish surveys of the Galactic plane have 
been carried out at 5 GHz with a resolution of a few ar- 
cminutes that is totally inadequate to complement the 
new generation of IR surveys (e.g., AltenhofT et al. 1978; 
Haynes et al. 1978; Gregory & Condon 1991; Griffith et 
al. 1994). Most higher resolution interferometric surveys 
have been carried out at the relatively low frequency of 
1.4 GHz (Zoonematkermani et al. 1990; Becker et al. 
1990; Condon et al. 1998; Giveon et al. 2005). Although 
the multi-configuration MAGPIS survey (Helfand et al. 
2006) and MGPS-2 survey (Murphy et al. 2007) are 
good for tracing the optically thin thermal emission from 
extended, evolved Hii regions and non-thermal super- 
nova remnants, they are not appropriate for the study 
of dense, thermal sources due to the unfavourable spec- 
tral index and compactness. Becker et al. (1994), with 
infills by Giveon et al. (2005) and White et al. (2005) 
have surveyed part of the inner galaxy (—10° < I < 42°, 
l&j < 0.4°) at 5 GHz with the VLA in C configuration giv- 
ing a resolution of 4"x9". However, these surveys cover 
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only one fifth of the northern GLIMPSE region and the 
spatial resolution is poorer than that of Spitzer IRAC, 
which causes problems with source identification. At this 
resolution the vast majority of their sources are unre- 
solved and thus we have no morphological information. 
In the southern hemisphere the AT20G survey at 20 GHz 
and '^30" resolution (Murphy et al. 2010) that was pri- 
marily designed to study the extragalactic radio source 
population is detecting many Hii regions, but again not 
at sufficient depth or resolution to address the science 
opened up by current surveys of the Galactic plane at 
infrared wavelengths. 

2. Scientific Motivation 

In this section we describe the various science topics 
that are addressed by the CORNISH survey. These cover 
the full range of galactic science from the formation and 
evolution of stars through to background extra-galactic 
sources. 

2.1. High-Mass Star Formation 

Our main science driver is the study of the formation of 
massive stars. These stars impart large amounts of UV 
radiation, wind energy, strong shocks and chemically en- 
riched material into the ISM during their lifetime. This 
is particularly true in starbursts when they form in dense 
concentrations (Leitherer 1999). However, relatively lit- 
tle is known about the physics that controls the forma- 



tion of stars more massive than about 8 Mq (10^ Lq) and 
hence the upper initial mass function (IMF). The fast col- 
lapse and contraction time scales, and extreme conditions 
resulting from increased turbulence and radiation pres- 
sure makes the theoretical treatment of their formation 
more challenging than for solar- type stars (e.g., Behrend 
& Maeder 2001; Yorke & Sonnhaiter 2002; McKee & Tan 
2003; Bonnell & Bate 2006). Their rarity and predilec- 
tion to form in dense clusters makes observational studies 
more challenging and high resolution a pre-requisite (e.g., 
Beuther et al. 2002). 

Current numerical simulations indicate that accretion 
via a disc can continue despite the strong radiation pres- 
sure (Krumholz et al. 2009; Kuiper et al. 2010) whilst 
sub-millimetre interferometry reveals evidence of flat- 
tened, disc-like structures (e.g., Patel et al. 2005; Tor- 
relles et al. 2007) around massive young stellar objects 
(MYSOs). Although these deeply embedded IR objects 
are already very luminous, they have yet to ionize the 
surrounding ISM, most likely due to ongoing accretion 
keeping them swollen and cool (Hoare & Franco 2007; 
Hosokawa & Omukai 2009; Davies et al. 2011). MYSOs 
drive highly coUimated ionized jets (Marti et al. 1998; 
Curiel et al. 2006) and equatorial disc winds (Hoare et 
al. 1994; Hoare 2006; Gibb & Hoare 2007) which are 
manifest as compact thermal radio emission (a few mJy 
for nearby examples) with a spectral index close to +0.6. 
What drives these jets and winds at speeds of a few hun- 
dreds of kms~^ (e.g., Bunn et al. 1995) is not known. 
A possible scenario is that the radio jets are accelerated 
and coUimated by an MHD mechanism as in low-mass 
YSOs early on, and these give way to less-coUimated ra- 
diation pressure dominated flows (Drew et al. 1998) at 
later stages, akin to the scheme outlined by Beuther & 
Shepherd (2005). 

The strong radio emission from H ii regions, formed 
once the star has contracted down to the main sequence 
and heated up, is a very useful probe of massive star for- 
mation (Hoare et al. 2007). It is unaffected by extinction 
and can thus be used to measure the formation rate of 
OB stars right across a galaxy (e.g., Condon, Cotton & 
Broderick 2002). The radio brightness of H ii regions can 
be used to infer the flux of Lyman continuum photons 
from the exciting star or cluster and hence constrain the 
spectral type and mass of the star(s). H ii regions provide 
one of the clearest signatures of spiral structure which 
can be related to the global distribution of molecular 
and atomic hydrogen from complementary surveys. This 
is particularly important in the context of the emerging 
picture of our Galaxy as having two stellar arms linked 
to the central bar whilst having four gaseous star forming 
arms (Benjamin 2009). 

High resolution studies of Galactic objects have re- 
vealed both large numbers of ultra-compact objects 
(UCHiis) and a preponderance of cometary morpholo- 
gies (Wood & Churchwell 1989; Kurtz et al. 1994; Walsh 
et al. 1998). These studies spawned a re-examination of 
the dynamical models of H ii regions (Van Buren & Mac 
Low 1992; De Free, Rodriguez & Goss 1995; WilUams, 



Dyson & Redman 1996). The latest cometary models 
(Arthur & Hoare 2006) require the young OB star to 
be located in a density gradient to explain the cometary 
shape, together with a strong stellar wind to give the 
limb-brightened appearance and an element of proper 
motion up the density gradient to reconcile the ionized 
and molecular velocity structure. Such a picture is con- 
sistent with one where the formation of these massive 
stars was triggered by the passage of an external shock- 
wave through the molecular cloud. 

There are many instances where the expansion of 
H II regions appears to be triggering formation of new 
episodes of massive star formation (e.g., Deharveng et 
al. 2005; Urquhart, Morgan & Thompson 2009; Moore 
et al. 2007; Zavagno et al. 2010; Thompson et al. 
2012). This is readily apparent in the 8/um images 
from the GLIMFSE survey, where the strong Polycyclic 
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Fig. 1. — Spectral energy distribution from the 
near-IR to the radio of an example UCHii region 
(G042.4343-00.2597, L = 1.8xlO*L0, d = 5.0kpc, solid 
circles) and MYSO (S140 IRS 1, L = 0.85x 10^ L©, 
d = 0.76 kpc, open circles). Sources of the IR data for 
the Hit region can be found in Mottram et al (2011) and 
the radio data come from Wood & Churchwell (1989) 
and Garwood et al. (1988). For S140 IRS 1 the IR data 
points are discussed in Maud et al. (2012) and the ra- 
dio points arc from Gibb & Hoare (1997) and Schwartz 
(1989). The flux densities from these objects have been 
reduced to place them at 15 kpc to reflect early B stars 
on the far side of the galaxy. The 2 mJy (5(t) sensitivity 
limit for the CORNISH survey is marked along with the 
5f7 limits for the UKIDSS, GLIMPSE, MIPSGAL, Her- 
schcl and JCMT surveys. Note that the sensitivities are 
well matched to the SEDs of our main targets and that 
the radio data are vital for distinguishing between the 
radio-loud UCH us and radio-quiet MYSOs. 



Aromatic Hydrocarbon (PAH) emission delineates the 
Photo-Dissociation Region (PDR) at the interface be- 
tween the ionized region and the molecular cloud. Fur- 
thermore the 24 iim images from MIPSGAL show up the 
warm dust heated by La within the ionized zone and cor- 
relates well with the radio free- free emission (e.g., Wat- 
son et al. 2008). Younger phases such as mid-IR bright 
MYSOs or maser sources are often seen just ahead of the 
ionization front sometimes within IR dark regions (Hoare 
et al. 2007). What fraction of massive star formation 
occurs in this way, as opposed to spontaneous collapse 
is unclear. Conversely, the action of outflows, expand- 
ing Hii regions and stellar winds can also disperse the 
remaining molecular cloud and quench further star for- 
mation (e.g.. Franco, Shore & Tenorio-Tagle 1994). How 
these opposing feedback mechanisms affect the resultant 
star formation rate and efficiency is unknown. 

Many aspects of massive star formation studies suffer 
from the lack of a large and unbiased sample of both 
UCH IIS and MYSOs. Samples based on the IRAS Point 
Source Catalogue (e.g., Kurtz et al. 1994; Molinari et al 
1996; Sridharan et al. 2002) are heavily biased towards 
bright, isolated sources due to the poor spatial resolu- 
tion of IRAS, which resulted in confusion in the densest 
regions of the plane, precisely where massive stars form. 
The situation for mid-IR bright MYSOs has been allevi- 
ated by the Red MSX Source (RMS) survey (Lumsden 
et al. 2002) that utilized the 18" resolution mid-IR sur- 
vey by the MSX satellite (Price et al. 2001) and high 
resolution, multi-wavelength follow-up to yield a large, 
well-selected sample of MYSOs. This is fairly complete 
out to about 15 kpc for O star luminosities (Urquhart et 
al 2011; 2012), and the GLIMPSE survey can easily push 
this down to the B star range (Robitaille et al. 2008). 
However, the IR spectral energy distributions of MYSOs 
and UCH us are very similar and therefore radio contin- 
uum data are vital to distinguish these two evolutionary 
phases (Figured]) (e.g., Urquhart et al. 2007). Hence, 
radio observations are crucial in distinguishing the 'radio- 
quiet' MYSOs from the 'radio-loud' UCHii regions. The 
large numbers of very red objects in the GLIMPSE sur- 
vey with colours of embedded YSOs makes radio follow- 
up of each one impractical. The uniform coverage of the 
CORNISH survey makes this automatic. 

The RMS survey also found many UCH ii regions, but 
such IR selected samples are biased against very em- 
bedded objects such as those found in IR dark clouds 
(IRDCs) (Egan et al. 1998). Many examples of methanol 
maser sources (EUingsen 2006), extended green objects 
(EGOs) (Cyganowski et al. 2008) and faint 24 ^m MIPS- 
GAL sources (Rathborne et al. 2010) are showing up 
in IRDCs. Is there radio continuum emission associ- 
ated with these very early phases of massive star for- 
mation? Cyganowski et al. (2011) performed deep ra- 
dio continuum observations of a sample of EGOs and 
found mostly weak or non-detected sources consistent 
with jets or winds in the diagnostic plots of Hoare & 
Franco (2007). The rarity of known MYSOs with lu- 
minosities above 10^ Lq raises the question of whether 



early O stars ever go through a detectable mid-IR bright 
MYSO phase or whether they manifest themselves firstly 
as UCH IIS. 

The clean selection function of the CORNISH survey 
combined with dynamical models and population simu- 
lation will allow us to determine the lifetime of this phase 
as a function of luminosity (e.g., Davies et al. 2011). To- 
gether with similar work for MYSOs from the RMS sur- 
vey and earlier phases from the Hi-GAL and SCUBA2 
surveys this will yield the first comprehensive picture of 
the evolutionary scheme from high mass star formation. 
Modelling the flux, size and shape distributions of the 
UCH II region population will test the dynamical models 
of UCH II regions and their relative locations in the dense 
molecular cloud material will enable statistical tests of 
triggering scenarios. These tests are only possible from 
uniform, wide-area surveys with sufficient sensitivity and 
resolution. 

In the high-mass star formation context, a spatial res- 
olution of about 1" is needed in a complementary radio 
survey for several reasons. Firstly, the IR source den- 
sity in regions of massive star formation is very high in 
dense star forming regions. Hence, to accurately identify 
the IR counterpart to a compact radio source, arcsecond 
resolution is needed. Secondly, UCH lis themselves often 
appear in complexes with separations of a few arcseconds 
(e.g., Sewilo et al. 2004). Finally, as previous surveys 
have shown, to reveal the structure of most UCH us, a 
resolution finer than an arcsecond is needed. 

2.2. Planetary and Proto-planetary Nebulae 

Zijlstra & Pottasch (1991) estimated that there should 
be about 23 000 PN in the galaxy, although only just 
over 3000 have been identified (Frew & Parker 2010). 
Combined radio and IRAS searches have been success- 
ful in turning up new PN (Van de Steene & Pottasch 
1995; Condon, Kaplan, Terzian 1999). A high resolution 
radio survey is particularly suited to discovering young, 
compact PN that are heavily reddened by line-of-sight 
extinction in the Galactic plane and we expect about a 
thousand of these in the survey area. Their radio flux 
densities are expected to be in the range of 5 — 50mJy 
for typical distances with sizes of a few arcseconds. A 
complete sample over a large region of the Galactic plane 
will constrain the density of PN and hence their forma- 
tion rate and distribution in the Galaxy. This has impli- 
cations for models of post-AGB evolution and mass-loss 
which are not well constrained (e.g., Casassus & Roche 
2001). 

As in the case of UCH us, the immediate progenitors 
of PN - the proto-PN (PPN) are also radio quiet since 
the central star has not yet reached 30 000 K. The radio- 
quiet PPN will also have very similar mid-IR colours to 
the radio-loud PN and so the only way to distinguish 
them is through their lack of radio emission. This is an- 
other illustration of the usefulness of non-detections as 
well as detections in the radio survey. PPN hold the key 
to understanding the envelope ejection process and test- 
ing the interacting winds models for the shaping of PNe 



(e.g., Mellema & Frank 1995). It is not known whether 
aspherical mass-loss at the end of the AGB is due to 
AGB star spin-up, detached or common-envelope bina- 
ries (Mastrodemos & Morris 1999) or dynamically impor- 
tant magnetic fields (Gardiner & Prank 2001). If we take 
the birth rate of PN to be about 1 per year (Kwok 2000) 
and the same for PPN with a PPN hfetime of 1500 years 
then we expect about 250 PPN in our survey region. 

2.3. Evolved OB stars 

Mass-loss from the later evolutionary stages of massive 
stars (eg. OB supergiants, Wolf-Rayet stars, and Lumi- 
nous Blue Variables) is of great importance in the evolu- 
tion of the ISM in our own and other galaxies, through 
both the kinetic power deposited by their stellar winds 
and their strong ionizing radiation fields. Thermal ra- 
dio emission plays a key role in determining mass-loss 
rates for the evolved massive stars since it arises at large 
distances from the underlying stars where the outflow 
has attained a constant velocity (e.g., Wright & Barlow, 
1975). Renewed interest in this important input to the 
ISM has arisen due to dumpiness of the winds (Fullerton 
et al. 2006). Radio deduced mass-loss rates may also pro- 
vide evidence of the bi-stability jump, a drop in ratio of 
wind-terminal velocity and the escape velocity (Lamers, 
Snow & Lindholm 1995), and perhaps important to un- 
derstanding of the LBV phase of stellar evolution (Vink 
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Fig. 2. — Illustration of the pointing pattern for a typical 
8 hr observing block. Each dot represents a pointing cen- 
tre. Scans in RA started at the bottom of the block and 
moved upwards before repeating to improve the uv cov- 
erage. The diagonal lines represent b = —1 and b = +1 
whilst the horizontal lines indicate the borders of this 
block. Note that pointings extended up to 6.4' outside 
of the nominal survey area to ensure uniform sensitivity 
within it. 



& de Koter 2002). 

Many of these massive stars exhibit non-thermal emis- 
sion in addition to the thermal emission from their stel- 
lar winds (Bieging et al. 1989), shown to be associated 
with massive binary systems, with the emission arising 
where the stellar winds of the two companions collides. 
The archetype of these colliding- wind binaries is WR140, 
where the non-thermal emission has led to the deter- 
mination of all the orbit parameters, the distance and 
the wind-momentum and mass ratios (e.g. Dougherty 
et al. 2005; Dougherty et al. 2011). Among WR stars 
that exhibit non-thermal emission, over 90 % are binaries 
(Dougherty & Williams, 2000), and there is now evidence 
that a similar fraction of non-thermal emitters is also ob- 
served in 0-star binaries (de Becker 2007). Challenges 
remain for the colliding-wind model. Not all stars with 
non-thermal emission are found to be binaries, and many 
binaries exhibit thermal emission, which among short- 
period systems might be attributed to absorption of the 
synchrotron power. Alternative mechanisms for acceler- 
ating electrons to necessary relativistic speeds have been 
proposed e.g. shocks within the stellar winds, that work 
in single star scenarios. 

Studies of both the thermal and non-thermal emitters 
are limited by the small sample size of the detected sys- 
tems, and the samples are often biased via spectral type, 
optical luminosity, or known binary systems. A large, 
unbiased sample of radio detected evolved OB stars has 
previously not been available. Since they inhabit the 
densest parts of the plane, there are likely many evolved 
OB stars hidden from optical view by high extinction. 
For a 5 GHz detection limit of 2 mJy, the limiting ra- 

















^^^^^^^^ '""'^^'^ 















310 300 290 280 270 260 

RA (dega) 



Fig. 3. — The entire survey area with the locations of 
the three phase calibrators indicated. The circles have 
radii of 10° showing that these calibrators were sufficient 
to cover the entire survey area. 



dio luminosity is 2.4 X 10^^(L'/kpc)^ ergs ^ Hz ^. Using 
the empirical relationship between Mand Lboi from Gar- 
many & Conti (1984) and the relationship between Lecm 
and Mfrom radio flux from Bieging et al. (1989) we 
expect this survey to detect all stars with a bolometric 
luminosity Lboi > 2 x W^Lq within 1 kpc. 

2.4. Active Stcirs 

Sensitive radio observations have revealed the ubiquity 
of high energy processes across the radio H-R diagram 
(Giidel 2002). As well as dynamo effects in convective en- 
velopes, high-resolution radio observations have revealed 
large-scale magnetospheric structures in various stellar 
types, including both low- and intermediate-mass pre- 
main sequence stars. Radio observations reveal dozens of 
such objects in Orion, and Chandra now finds over 1000 
active stars in Orion (Garmire et al 2000). With 5 GHz 
flux densities reaching several tens of mJy, these sources 
will be detected to several kpc. Particularly intriguing is 
the role of magnetic fields in intermediate-mass pre-main 
sequence stars, which are not expected to have convective 
envelopes and hence large scale dynamos (e.g., Garring- 
ton et al 2002). Combined radio, IR and X-ray studies 
will reveal how common these magnetic phenomena are. 

2.5. Active Binetries 

Mass transfer in close binaries often leads to strong 
activity in radio and other wave-bands. RS CVn stars, 
cataclysmic variables and others are likely to be de- 
tected in the survey. However, the most important 
class of active binary are the so-called micro-quasars 
where the accreting object is a black-hole or neutron 
star (Mirabel & Rodriguez 1999; Fender 2002). These 
produce relativistic jets analogous to AGN where ejec- 
tions of plasma can be followed in real time (Fender et 
al. 1999). These objects are perhaps the best laborato- 
ries for studying the relativistic jet formation process and 
test magneto- hydrodynamic models (e.g., Meier, Koide 
& Uehida 2001). The key signature of these jets is strong 
synchrotron emission in the radio band and, since they 
often originate from massive stars, one expects to find 
them close to the Galactic plane and in the locale of gi- 
ant molecular clouds. In such dense environments the 
X-ray and optical signatures of accretion may be heavily 
absorbed, preventing the detection of accretion by more 
conventional means. Hence, a co-ordinated radio and IR 
survey should uncover many new examples of this phe- 
nomena. 

There is a strong relationship between the X-ray state 
and radio behaviour in these objects, which has led to 
the concept of disk-jet coupling. Several black hole can- 
didate objects have flat or inverted radio spectra when 
in the low/hard X-ray state, indicative of the presence 
of a self-absorbed compact jet supported by a hot X- 
ray corona. The jets have been resolved by VLBI in 
GRS1915-I-105 (Dhawan et al. 2000) and Cygnus X-1 
(Stirling et al 2001). When the accretion rate increases 
the thermal disk becomes dominant and the X-rays be- 
come softer. Radio outbursts often occur during the 



transition between hard and soft states (Fender, Homan 
and Belloni 2009) due to the formation of a strong shock 
in the out-flowing jet, and radio interferometer observa- 
tions are able to follow the evolution of the flow (e.g.. 
Fender et al 1999). Our survey has the resolution to 
resolve such structures and at high accretion rates the 
objects are detectable throughout the disk of the galaxy. 

2.6. Extragalactic Sources 

There will be a significant number of background ex- 
tragalactic sources detected during the survey. Using 
equation A2 in Anglada et al. (1998), which is based 
on the source count data from Condon (1984), at our 
completeness level we expect to detect about 2500 extra- 
galactic sources in the 110 square degree area. The vast 
majority of these radio galaxies will be at too high a red- 
shift to determine distances via H I detections with cur- 
rent telescopes (Magliocchctti et al. 2000; Obreschkow 
et al. 2009). As well as penetrating the deepest part of 
the zone of avoidance these sources will be useful in other 
respects. The most compact ones are potential phase cal- 
ibrators for use in future interferometric studies in the 
Galactic plane. Being at 5 GHz the CORNISH survey 
will preferentially pick up flat spectrum sources, which 
are potential mm interferometric phase calibrators. Ex- 
tragalactic sources that are also visible in the optical or 
near-IR can be used to bring together the radio and op- 
tical coordinate reference frames which is important for 
multi-wavelength high resolution studies in the Galactic 
plane. 

Radio surveys can also deliver a set of extragalactic 
background sources that can act as pencil beam probes 
of molecular clouds in follow-up absorption line studies 
(Greaves & Nyman 1996). Sufficiently strong sources 
could be used for Zeeman absorption line measurements 
of the magnetic field strength (Crutcher 1999). For 
instance, radio sources that lie behind dense, starless 
molecular cloud cores can be used to probe the initial 
conditions for star formation in these dense clouds with 
future sensitive instruments such as the SKA. Polarized 
background sources can also be used to probe the mag- 
netic field of the Galaxy, individual H ii regions and su- 
pernova remnants via Faraday rotation (Taylor et al. 
2009). High frequency data are useful in disentangling 
components in sources with complex Faraday depth spec- 
tra (O'SuUivan et al. 2012). 

3. Survey Specification and Implementation 

The scientific motivation described in the preceding 
section required a radio continuum survey at around 
1 mJy sensitivity, 1" resolution and at a frequency sig- 
nificantly higher than the 1.4 GHz traditionally used for 
surveys. Given the positive spectral index of many of the 
target sources then the higher the frequency the better. 
For VLA receivers the 8 GHz band would normally be the 
optimum choice in terms of sensitivity for objects with 
spectral index around -1-0.6. However, due to the pri- 
mary beam size, the time required to survey a given area 
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Fig. 4. — Example data from the CORNISH survey alongside three-colom images from the Spitzer GLIMPSE (3.6 /im, 
4.5 ^m and 8.0 //m) and UKIDSS (JHK) surveys. Top: A typical cometary Hii region, G032. 1502+00. 1329. Note the 
good correspondance between the radio emission and the brightest 8 fim emission in the GLIMPSE image. There is 
also strong extinction ahead of the cometary as expected if the OB star formed in a density gradient. Middle: A 
candidate PN that was previously unknown, G012. 3830+00. 7990. It is clearly seen as an isolated source at 8 fim with 
a very faint counterpart in the near-IR. In the longer wavelength Spitzer MIPSGAL data (not shown here) it is bright 
at 24 yum but fainter at 70 /im unlike H ii regions that are brighter at 70 /im. Bottom: A previously unknown radio 
star, G030. 2357— 00.5719, that has blue colours in the near- and mid-IR. 



increases as ly^. Thus a large survey is much more prac- 
tical at 5 GHz. The VLA 5 GHz sensitivity is similar to 
that at 8 GHz, but the survey only takes half the time. 
For a frequency of 5 GHz the VLA in B configuration 
delivers a resolution of 1.5" that satisfies the scientific 
requirement. A significant fraction of the area to be sur- 
veyed is below declination —15° where it is necessary to 
switch to the BnA configuration to ensure that the beam 
elongation does not degrade the resolution. As with any 
single configuration snapshot survey the spatial dynamic 
range is limited and structures larger than about 12" in 



size are not imaged well in the CORNISH survey. 

In order to image the full (8.9') primary beam with- 
out image degradation due to bandwidth smearing these 
observations were carried out in spectral line mode. A 
25 MHz bandwidth in 4IF mode was used to maximise 
the sensitivity. This has eight 3.1MHz channels which 
degrades the peak response by only a few percent at the 
edge of the primary beam. 

To achieve a noise level of 0.3mJybeam~^ we required 
a total on-source integration time of 90 seconds per point- 
ing. This noise level is similar to that used by Wood & 
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Fig. 5. — Further example data. Top: A previously unknown radio star, G029. 7188— 00.0316, that has red 
colours in the near- and mid-IR. Middle: A candidate radio galaxy where we only see one unresolved component, 
G049. 6617— 00.0543. No counterpart is seen in the IR. Bottom: A candidate radio galaxy where we see two unresolved 
components, G026. 7179— 00.8290. Again, no counterpart is seen in the IR. 



Churchwell (1989) and Kurtz et al. (1994) to detect and 
map UCHii regions throughout the Galaxy. If we use 
equations 1 and 3 from Kurtz et al. for UGH lis with an 
unresolved flux density of 2 mJy placed 20 kpc away we 
find that we can detect all objects powered by a B0.5V 
star or earlier. 

To achieve a sensitivity that is uniform across the sur- 
vey to within 10 %, we adopted the NVSS hexagonal pat- 
tern of pointing centres (Gondon et al 1998). Scaling the 
NVSS pointing pattern to 5 GHz gives a hexagonal grid 
separation A of 7.4'. The observations were carried out 
in the pseudo-holography mode developed for the NVSS 
survey (see Gondon et al. 1998), scanning in RA be- 
tween pointings. The separation between pointings in 
RA was 7.4' = (2 sin 30°) A whilst those in Dec were sep- 



arated by 6.4'= (sin60°)A. The pointing pattern for a 
typical observing block is shown in Figure [21 Forty such 
blocks were completed in 360 hours of time to complete 
the survey. The survey was spread over two trimesters, 
with 21° < ^ < 49° being observed in the third quarter of 
2006 and the rest in the last quarter of 2007 and begining 
of 2008. 

Three phase calibrators (1832-105, 1856+061 and 
1925-1-211) were used to cover the whole survey with 
all blocks being within about 10° of one of these (see 
Figure [3|) . This continuity of phase calibrator allowed 
additional cross checks of the overall calibration of the 
survey. A primary fiux calibrator (1331-1-305) was ob- 
served at the beginning and a secondary flux calibrator 
(0137-f 331) at the end of each 8 hour observing nm. For 
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Fig. 6. — Plot showing the dweh time achieved and the final area covered by the CORNISH survey. 
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Fig. 7. — Plot showing the efi^ective number of antennas achieved during the CORNISH survey after flagging. 



the BnA observations 9 rows were completed in shorter 
6 hour runs. 

Complete coverage of the survey area was achieved and 
indeed extended to / = 65.4° that was the eventual limit 
of the GLIMPSE survey. A log of the observation dates 
can be found in Table [5] Example images, representative 
of the different types of science discussed in section [2l are 
presented in Figures |4] and El alongside those from the 
GLIMPSE survey. The complementarity of the similar 
resolution radio and IR datasets is very apparent. 

Figure [H] shows the dwell time across the survey area 
where the uniform coverage is clearly seen. About 5 % 
of the survey pointings were repeated due to weather or 



technical problems. Only one scan had a single 45 second 
pass rather than the usual two (see below) because the 
repeat scan also failed. Along the edges at 6 = — 1° and 
l°on the end points of scans the on-source time varied 
due to the nature of the raster scanning technique. Dur- 
ing the survey the VLA antennas were being retro-fitted 
into EVLA antennas, which meant that there were often 
significantly less than the nominal 27 antennas. In addi- 
tion to the usual flagging of antennas and baselines due 
to technical issues this meant that the effective number of 
antennas varied across the survey as shown in Figure [71 
The median noise level achieved over the survey area 
was 0.33mJy beam^^ (Purcell et al. 2012) very close to 
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Fig. 8. — Top left: Simulated VLA beam for a single 2-mimitc snapshot at an hour angle of —4 hours and declination 
—3°. Top right: Simulated VLA beam for two 1-minute snapshots at hour angles of —4 hours and —2 hours, i.e. 
separated by 2 hours. Bottom left: Simulated VLA beam for two 1 minute snapshots at hour angles of —4 hours 
and hours, i.e., separated by 4 hom-s. Bottom right: Simulated VLA beam for two 1 minute snapshots at hour 
angles of —4 hours and +2 hours, i.e., separated by 6 hours. 



Table 2 
Log of observations. 



Dates 


/ ranges observed 


12/07/2006 to 16/09/2006 


21.1° 


to 48.9° 


28/09/2007 to 06/10/2007 


10.0° 


to 16.1° 


27/10/2007 to 04/02/2008 


16.1° to 21.1° 


and 48.9° to 65.4° 
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Fig. 9. — The variation in the synthesised beam major axis, minor axis and position angle over two adjacent observing 
blocks covering 21.5° < / < 24°. 
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Fig. 10. — Plot showing the major beam axis achieved throughout the survey. 
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that expected. To determine the nature of the sources 
in the CORNISH survey it will be very useful to evalu- 
ate the spectral index, which requires data at other radio 
wavelengths. The principal existing source of these data 
will be the 1.4 GHz plane surveys. For the survey car- 
ried out using the VLA in B configuration by White et al. 
(2005) that overlaps with 91 % of the CORNISH survey 
area their quoted median noise level is 0.90 mJy beam" ^, 
although there is significant variation over the area. This 
means that typically a source would need to have a spec- 
tral index steeper than —0.8 to be guaranteed to be de- 
tected in both. If we take the quoted completeness limit 
of 13.8 mJy from White et al. (2005) then optically thin 
Hii regions with a = —0.1 will need to be stronger than 
12 mJy in CORNISH to be also detected in the 1.4 GHz 
survey. For optically thick Hii regions with a = 2.0 this 
rises to 150 mJy. For 55% of the CORNISH area muhi- 
configuration VLA 1.4 GHz data exists from the MAG- 
PIS survey (Hclfand et al. 2006) with a noise level that is 
more uniform and three times lower than the White et al. 
(2005) survey. Therefore, for most non-thermal sources 



in the CORNISH survey with a typical a = —0.6 then 
sufficient data exists to characterise their spectral index. 
However, for thermal sources with a positive spectral in- 
dex follow-up observations will be required to measure 
the spectral index. 

4. Beam Shape 

To improve the uv coverage and reduce the sidelobes 
the 90 second integrations were split into two separate 
snapshots of 45 seconds, each separated by 4 hours in 
hour angle. This strategy was arrived at after simu- 
lations were made of the beam shape for various hour 
angle separations. Figure [8] shows the simulated beams 
for a single 2 minute snapshot as well as two one-minute 
snapshots separated by 2, 4 and 6 hours. Table|3]lists the 
beam parameters for each of these simulations. The two 
snapshots clearly give better beam structure, with the 
ones separated by 6 hours giving the best beam shape. 
However, we used a 4 hour separation for the survey as 
this made implementation and scheduling easier in the 8 



Table 3 

Simulated VLA beam parameters for the different snapshot strategies shown in Figure [H 



HAi (hrs) 


HA2 (hrs) 


femaj (") 


femin (") 




-4 




3.3 


1.4 


2.3 


-4 


-2 


2.5 


1.4 


1.8 


-4 





2.3 


1.4 


1.6 


-4 


+2 


2.0 


1.5 


1.3 



1.5F r 




Galactic Longitude (deg) 



Fig. 11. — Plot showing the minor beam axis achieved throughout the survey. 
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Fig. 12. — Plot showing the beam position angle achieved throughout the survey. 



hour observing blocks when the Galactic plane was ob- 
servable. 

A typical scanning row contained 20 pointings (see Fig- 
ure [2]), which including the drive time of about 5 seconds 
between pointings took about 17 minutes. Following each 
scan row a phase calibration of 2 minutes on-source was 
carried out. Eleven of these rows formed an observing 
block of about 3.6 hours, after which the whole block 
was repeated, giving the second snapshot at the appro- 
priate 4 hour separation in hour angle. For the BnA 
observations the observing runs were only 6 hours long 
due to the lower elevation and the two snapshots were 
then only 3 hours apart. 

The actual beam shape achieved in two consecutive 
blocks from survey observations carried out during 2006 
is shown in Figure El This is the synthesised beam (im- 
aged with a robust parameter of zero) after the combi- 
nation of the two individual 45 second snapshots. The 
major axis of the beam is slightly worse at the beginning 
of the blocks (low declination end) than that predicted in 
Table[3]for a similar declination. This is because schedul- 
ing constraints meant that the 2006 blocks did not quite 
transit mid-way through the 8 hour observing session 
which resulted in a more asymmetric beam than other- 
wise achievable. The beam at the end of the block (high 
declination end) should also have a major axis of about 
2.3" from a combination of the first snapshot at hour 
angle and a second one at hour angle -1-4. Another 
manifestation of the asymmetry in the scheduling is in 
the beam position angle that rotated from —40° to -1-20° 
throughout the block. The minor axis of the beam varied 
less across the blocks, as expected from Table [3l 

The beam parameters over the full survey area are 
shown in Figures [10] to [121 Observations conducted dur- 
ing the 2007/2008 were scheduled more symmetrically 
than in 2006, although the inevitable variation in beam 
shape is still present. The large change around / = 16° 
and below is due to the switchover to the BnA config- 
uration for declinations below —15°. Here the position 
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Fig. 13. — Distribution of major (grey) and minor (left- 
slanted) beam axes. 



angle of the elongated beam changes from being mostly 
north-south to an east-west direction. 

Histograms of the beam parameters achieved are shown 
in Figures [T3] and [T31 The median major axis size 
achieved was 1.5" with a standard deviation of 0.3" and 
a tail up to about 3". As expected the minor axis size 
has a fairly narrow distribution around 1.2" apart from 
the small number of higher resolution observations due to 
the BnA configuration data. The distribution of beam 
elongations shows that 98 % of fields have elongations 
less than 2 and 74% less than 1.5 in agreement with the 
predictions of Table [31 This shows that in the main the 
adopted strategy of two snapshots separated by 4 hours 
in hour angle did deliver the desired constraints on beam 
shape. 

5. Conclusions 

The CORNISH survey is a new sensitive, high spatial 
resolution Galactic plane survey at 5 GHz. It is targeted 
at compact thermal sources and UCH ii regions in par- 
ticular. Conceived to complement the Spitzer GLIMPSE 
survey, it will systematically address key questions in 
massive star formation such as the lifetime and evolution 
of the UCH II phase as a function of exciting star parame- 
ters and environment. Uniform coverage of the northern 
GLIMPSE region enables the distinction between radio- 
loud and radio-quiet objects such as MYSOs/UCHiis 
and PPN/PN that otherwise have very similar mid-IR 
colours. Legacy survey science in combination with other 
recent and upcoming optical, near-IR, far-IR, sub-mm 
and longer wavelength radio surveys of the plane will be 
possible with applications in evolved stars, active stars 
and active binaries. 

Data taking for the CORNISH survey has been com- 
pleted and the basic goal of uniformity of coverage 
achieved. The data reduction, source extraction and sta- 
tistical properties of the survey and sources will be dis- 
cussed in the forthcoming paper by Purcell et al. (2012). 
At that time the catalogue and image data will be pub- 



hcally available at http://www.ast.leeds.ac.uk/cornish 
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Fig. 14. — Distribution of beam elongation. 



We have also achieved the aim of controlhng the beam 
shape across the survey such that the elongation of the 
synthesiscd beam is mostly less than 1.5 and nearly al- 
ways less than 2. This was achieved by using a strategy 
of two snapshots separated by 4 hours in hour angle. A 
more uniform and rounder beam could be achieved by 
more snapshots but that would increase the overheads 
and has diminishing returns with the VLA. Our simula- 
tions show that snapshots separated by 6 hours in hour 
angle would have produced a better beam. However, 
this would have complicated the scheduling of the sur- 
vey into 8 hour runs and more than likely compromised 
the very uniform coverage that we ultimately achieved. 
Such tradeoffs will be the case for any survey with the 
EVLA at these declinations and for future Galactic plane 
surveys in particular. 

The National Radio Astronomy Observatory is a facil- 
ity of the National Science Foundation operated under 
cooperative agreement by Associated Universities, Inc. 
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